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Abstract 

7-ray bursts have baffled theorists ever since their accidental discovery at the 
sixties. We suggest that these bursts originate in merger of neutron star bina- 
ries, taking place at cosmological distances. These mergers release ~ lO^^ergs, 
in what are possibly the strongest explosions in the Universe. If even a small 
fraction of this energy is channeled to an electromagnetic signal it will be de- 
tected as a grbs. We examine the virtues and limitations of this model and 
compare it with the recent Compton 7-ray observatory results. 



Prologue: 7-Ray Bursts circa 1973 

7-ray bursts (grbs) were accidentally dis- 
covered ahead of their time. Had it not been 
for the need to verify the outer space treaty 
of 1967 (which forbade nuclear experiments 
in space) we would not have known about 
these bursts until well into the next century. 
No one would have proposed a satellite to 
look for such bursts, and had such a pro- 
posal been made it would have surely turned 
down as too speculative. The VELA satel- 
lites with omni-directional detectors sensitive 
to 7-ray pulses, which would have been emit- 
ted by a nuclear explosion, were launched 
in the mid sixties to verify the outer space 
treaty. These satellites never detected any nu- 
clear explosion. However, as soon as the first 
satellite was launched it begun to detect puz- 
zling, perplexing and above all entirely unex- 
pected bursts. The lag between the arrival 
time of the pulses to different satellites gave 
a directional information and indicated that 
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the sources are outside the solar system. Still, 
the bursts were kept secret for several years, 
until Kelbsdal Strong and Olson described 
them in a seminal paper|Q in 1973. 

The unexpected discovery sent theorist in 
a search for a source model. Specific models 
ranged from comets to cosmic strings. 

False Clues? 

The rapid fluctuation in the signal (less 
than 10ms) suggested a compact source, a neu- 
tron star or a black hole. Several other clues 
focused the attention of theorists towards neu- 
tron stars at the disk of the galaxy. 

First, came an analytic estimate p| of the 
optical depth to 77 — > e~^e~ . For an impulsive 
source we have: t^j ~ (TtPD"^ / R^meC^ where 
F is the fluence (~ lQ~^ergs/ crn^ in the early 
detectors and ~ lQ~^ergs/cm? in Compton- 
GRO), D is the distance to the source and R 
its size (expected to be less than 10^ from tim- 
ing arguments). Since r^^ > 1 for D > lOOpc, 
it was argued that the sources must be at 



the disk of the galaxy. Otherwise, it was ar- 
gued an optical thick system will cool down 
and radiate its energy in the x-ray uv or opti- 
cal band and not as 7-rays. The non-thermal 
spectrum also indicated that the sources are 
optically thin. Incidentally,, it was the con- 
frontation between this argument and the indi- 
cations from Compton-GRO that grbs are cos- 
mological (which we discuss later) that have 
lead to claims that grbs require "new physics" . 
We will see that it ain't necessarily so. 

A very strong and long (1000 sec) burst 
was observed on March 5th 1979. The posi- 
tion of the burst coincided with a SNR rem- 
nant in the LMC supporting the idea that grbs 
originates on neutron stars. 

Another clue came from the observation of 
absorption lines p, ^. The lines were inter- 
preted as cyclotron lines in a lO^^G magnetic 
field, a field strength that is found only on 
neutron stars. 

These clues and others have led to the con- 
sensus that grbs arise on neutron stars in the 
disk of the galaxy 0, possibly in their magne- 
tosphere. 

Bursts distribution circa 1991 

There were, however some indications that 
the sources might not be galactic. In 1975 
Usov and Chibisov suggested to use a 
logN-LogS test to check if the bursts have 
a cosmological origin. Later, in 1983 van 
den Bergh |^ analyzed the distribution of the 
46 bursts that were known at that time and 
from the isotropy of this distribution he con- 
cluded that the sources are either local at dis- 
tances of less than half of the galactic disk 
scale height or cosmological at redshift z > 
0.1 (See also |^). The cosmological solution 
was accepted with skepticism since with typ- 
ical fluencies of 10~^ergs/cm'^ the bursts re- 
quire lO'^^ergs if they originate at distances 
larger than lOOMpc! In 1986 Paczyhski M 



argued that the bursts are cosmological and 
suggested that some of the burst are lensed 
by intervening galaxies and that this will pro- 
vide an observational test to the cosmologi- 
cal hypothesis. In 1989 Eichler, Livio, Pi- 
ran and Schramm |IU| (see also [|TT|) suggested 



that the bursts originate in neutron star merg- 
ers at cosmological distances (The possibility 
that grbs might be produced in neutron star 
mergers was also mentioned without specify- 
ing a model by [0, |, 0, Q). However, in 
1991, just before the Compton-GRO results 
were announced, Atteia etal. [0 reported (at 
a 3(7 level) that the 244 bursts observed by 
the spacecrafts Venera 13, 14 and Phebus are 
concentrated towards the galactic plane, sug- 
gesting a disk population after all. 

Neutron Star Binaries 

Another seemingly unrelated and unex- 
pected discovery was make in 1975 by Hulse 
and Taylor 
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who found a pulsar, PSR 
1913+16, that was orbiting around another 
neutron star. No one has predicted that such 
systems exist, but in retrospect it was not sur- 
prising. More than half of the stars are in bi- 
nary systems. If some of these binaries survive 
the two core collapses (and the supernovae ex- 
plosions) needed to produce the neutron stars 
they will end in a binary pulsar. 

The binary pulsar have proven to be an 
excellent laboratory for testing General Rela- 
tivity. The binary system emits gravitational 
radiation which is too week to be detected di- 
rectly, but its backreaction could be observed. 
By carefully following the arrival time of the 
pulsar's signals Taylor and his collaborator 
have measured the pulsar's orbit. They have 
shown that the binary spirals in just in the 
right rate to compensate for the energy loss 
by gravitational radiation emission (with two 
neutron stars, tidal interactions and other en- 
ergy losses are negligible). For PSR 1913+16 



the spirahng in takes place on a time scale of 
tgr = 3 X lO^years, in excellent agreement 
with the general relativistic prediction |T7|]. 
These observations not only confirm the gen- 
eral relativistic prediction, they also assure us 
that the orbit of the binary is indeed decreas- 
ing and that inevitably in 3 x 10^ years the 
two neutron stars will collide and merge! 

Source Count and Event Rate 

For many years only one binary pulsar was 
known. A simple estimate based on the ob- 
servation of one binary pulsar in several hun- 
dred observed pulsars led Clark etal. 
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conclude that about 1 in 300 pulsars is in a 
binary. With a pulsars' birth rate of one in 
fifty years this led to a binary birth rate of 
one in 10"^ years. Assuming a steady state, 
this is also the merger rate. This estimate ig- 
nored, however, selection effects in the detec- 
tion of binary pulsars vs. regular ones. Specif- 
ically PSR1913+16 is an extremely bright pul- 
sar which is detectable from much larger dis- 
tance than an average pulsar. Currently there 
are four known binary pulsars and an analy- 
sis based on their luminosities and life times 
T9| , pO(] suggests that there are ~ 10"^ — 10^ 



neutron star binaries in the galaxy and that 
their merger rate is one per lO^years per 
galaxy. This corresponds to ~ 100 mergers 
per year in galaxies out to a distance of 1 
Gpc and about 10^ per year to the horizon. 



Narayan, Piran and Shemi |^ also predict 
that a similar or somewhat smaller population 
of neutron-star black hole binaries will exist. 

Neutron Star Mergers 

It was immediately realized, after the dis- 
covery of PSR1913-I-16 that the binary pro- 
duces a unique chirping gravitational radiation 
signal during the last seconds before the neu- 
tron stars merge. These signals are probably 



the best candidates for detection of gravita- 
tional radiation. However, these events are 
rare and to observe them (gravitationally) in 
our life time we must turn to extragalactic 
events. This is the aim of the advanced gravi- 
tational radiation detectors like LIGO pT|. 

As the strongest sources of gravitational 
radiation neutron star mergers attracted the 
attention of relativists, but most astronomers 
ignored them as being too rare to be of inter- 
est. Clark and Eardley have shown that 
the binding energy released in a neutron star 
binary merger is ~ 5x 10^^ — lO^'^er^'s, making 
these events possibly the most powerful explo- 
sions in the Universe. A significant fraction of 
this energy is emitted as gravitational radia- 
tion, both prior and during the collision. A 
very sophisticated gravitational radiation de- 
tector, LIGO, is built to detect these gravita- 
tional radiation signals. But it will be around 
the turn of the century when it is operational. 

As the neutron stars collide a shock forms 
and the stars heat up. Most of the binding en- 
ergy is emitted as neutrinos |^|. The neutrino 
burst is comparable or slightly stronger than a 
supernova neutrino burst (such as the one de- 
tected by Kamiokande and 1MB from 1987A). 
To detect extragalactic events at cosmological 
distances we need a detector which is ~ 10* 
times larger than those detectors. With regu- 
lar supernova neutrino bursts being a hundred 
times more frequent it is clear that these neu- 
trino signals are not the prime candidates for 
detection. 

Neutron star mergers are hiding from us 
by emitting their energy in two channels with 
extremely small cross sections. If even a small 
fraction of the energy is channeled to an elec- 
tromagnetic signal, its much large cross sec- 
tion will make it much easier to observe. For 
many years, I kept wondering what are the 
possible observational consequences of such 
events . 



Energy Conversion 
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Goodman, Dar and Nussinov 
gested that the neutrino-anti neutrino anni- 
hilation u + u — > e'^ + e~ converts a small 
fraction of the neutrino supernova burst to 
electron-positron pairs which in turn annihi- 
late to 7-rays, heat the surrounding envelop 
and provide the energy required to power the 
supernova shock wave. In 1989, Eichler, Livio, 
Piran and Schramm [1^ (see also [11|) sug- 
gested that the same mechanism operates in 
neutron star mergers and converts ~ 10"'^ of 
the emitted energy to pairs and 7-rays. This 
corresponds to lO^^ergs, roughly sufficient for 
detection of the bursts from cosmological dis- 



tances. Eichler etal. ITOi used the old estimate 
of Clark etal. 



||18| for the merger rate and sug- 
gested that these events would be detected by 
Compton observatory as grbs. 

More recently, alternative energy genera- 
tion mechanism such as magnetic field recom- 
bination 12311 or accretion onto the neutron star 



have been proposed and it was argued that 
they provide comparable amounts of energy. 

Fireballs and Relativistic Effects 

If grbs are indeed cosmological they are ini- 
tially optically thick, as Schmidt @ have ar- 
gued. How can there be a 7-ray burst from 
such a source? Goodman []T3| considered a 
dense sphere of 7-ray photons and pairs, which 
he called a fireball. He has shown that as long 
as the fireball is optically thick the radiation- 
pair plasma will behave like a fiuid with p = 
p/3. The fireball will expand and cool, just 
like the early Universe (unlike our Universe 
the gravitational force is unimportant). As 
the fireball cools its temperature drops with 
T (X 1/R until the electron positions annihi- 
late (the annihilation is complete at T ^ 20 
keV) and the radiation escapes. The radiation 
fiuid has reached in the meantime a relativistic 



velocity relative to an observer at infinity and 
its Lorentz factor F ^ Resc/Ro ~ T^/Tesc ~ 
10^ — 10^. The escaping photons, which have a 
typical energy of 20 keV in the local frame are 
blue shifted relative to an observer at infinity 
and their observed energy is eobs ~ TTe^c ~ Tq, 
of the same order as the initial energy. In this 
way the optical depth argument which limited 
the distances to the sources is bypassed and 
there is no need to introduce "new physics" to 
explain grbs from cosmological distances. 

Paczyhski ||^ have shown that similar ef- 
fects take place if the radiation is released in a 
quasi- stationary manner. In this case the ra- 
diation fiows out relativistic wind, with 
T (X 1/R and F cx i?. The radiation ceases to 
behave like a fiuid and escapes when T 20 
keV in the local frame. The escaping x-ray 
photons are blueshifted to much higher ener- 
gies in the observer frame. 

Do Fireballs Work? 

The fireball model faces two serious objec- 
tions: the origin of the observed nonthermal 
spectrum and the effects of baryons. 

There is no clear way to explain the non- 
thermal spectrum from a fireball that passes 
an optically thick phase and termalizes. It is 
possible that different regions in a realistic, in- 
homogeneous fireball move with significantly 
different Lorentz F factors and that the ob- 
served spectrum is a blending of thermal spec- 
tra to a non thermal one. Simple calculations 
of the spectrum of a spherical fireball |^ show 
some deviation from a thermal spectrum, but 
it is not large enough. Alternatively one could 
hope that the spectrum would become non- 
thermal in the transition from optically thick 
to optically thin regimes. However, this transi- 
tion takes place at ~ 20keV in the local frame. 
The energy injected from annihilation at this 
stage is insignificant and the temperature is 
too low for inverse Compton scattering to be 



effective ||2^. It seems that there is no clear 
mechanism that will modify the photons' black 
body spectrum in this stage. 

One expects that some baryons will be 
ejected into the fireball. Shemi and Piran P? 



tors 10 < r ^ Eq/Mc^ < 10^ Recently 



have shown that the baryons have two effects. 
For lO-i^M© < M < 10^^ Mq{Eo/ 10^^ ergs) 
the baryons dominate the opacity (long af- 
ter all the pairs have annihilated) without in- 
fluencing the fireball's inertia. The fireball 
continues to be optically thick until Tg = 
aTM/E? = 1. This leads to a longer accelera- 
tion phase and to a larger final Lorentz factor 
Tf ^ R/Rq ~ Tq/T. However, the final energy 
of the escaping radiation remains unchanged 
with e^TT ^ Tq. 

Larger baryonic load changes the dynamics 
of the fireball. As the fireball expands p oc R~^ 
while e oc r'^ If M > IQ-^ Mq{Eo/ 10^^ ergs) 
the baryonic rest mass will dominate the en- 
ergy density and the fireball's inertia before 
the fireball becomes optically thin. In these 
cases all the energy will be used to acceler- 
ate the baryons with Ek = Mc^T [Eq + 
Mc^)/{EoT/Mc^Tq + 1). The final outcome of 
a loaded fireball will be relativistic expanding 
baryons with F k, Eq/Mc^ and no radiation at 
all. 

Several ideas have been proposed to avoid 
the baryonic load problem. These include: (i) 
Separation of the radiation and the baryons 
due to deviations from spherical symmetry - 
the radiation escaping along the axis and the 
baryons being ejected preferably in the equa- 
torial plane \T^\ and (ii) generation of a radia- 
tion fireball with very small amounts of matter 



via magnetic processes 23 . 



Energy Conversion, Once More 

If the baryonic contamination is in the 
range IQ'^Eq < Mc^ < O.l^o all the initial 
fireball energy will be converted to extremely 
relativistic protons moving at a Lorentz fac- 



Meszaros, and Rees ||29| (see also [^) sug- 
gested that this energy could be converted 
back to 7-rays when this baryons interact with 
the surrounding interstellar matter. A shock, 
quite similar to a SNR shock, forms and it 
cools predominantly via synchrotron emission 
in the x-rays. The x-ray photons will be 
blueshifted to 7-rays in the observer frame due 
to the relativistic velocity of the fireball. The 
relativistic motion will also lead to a short time 
scale for the burst. Alternatively, the acceler- 
ated baryons could interact with a pre-merger 
wind that surrounds the fireball E^, |2H|. In 



both cases the interaction with the surround- 
ing material will lead once more to the conver- 
sion of the energy: from kinetic energy back to 
radiation. Since this phase is taking place in 
an optically thin region the photons will not 
thermalize and the emerging spectra will be 
non thermal, as observed. Thus, this process 
seems to resolve at one stroke both major ob- 
jections to the fireball scenario at one stroke. 

7-ray bursts distribution circa 1992 

The Compton 7-ray observatory was 
launched in the spring of 1991 (see ||30| for 
a review). It includes an omnidirectional 7- 
ray burst detector (BATSE) which, with a 
limiting sensitivity of ~ ergs/ am?, is the 
most sensitive detector of this kind flown. By 
the summer of 1992 BATSE has detected more 
than 400 bursts, more than all previous de- 
tectors combined. BATSE is also capable 
of obtaining a directional information on the 
bursts on its own. Within four month from its 
launch BATSE has collected enough data to 
conclude that the distribution of grbs sources 
is isotropic When the V/Vmax test was 

applied to the burst intensity distribution is 
was shown that the sources are not distributed 
homogeneously in space and that there is a 



concentration of sources towards us [31 . 



These two observations rule out all lo- 
cal galactic disk models. The observations 
are consistent with three possible populations: 
(i) Cosmological population (ii) Galactic halo 
population with a large core radius (> 50kpc) 
and (iii) A population, such as comets at the 
Oort cloud, centered around the solar system. 

We will turn to the second and third pos- 
sibilities, before summing up the status of the 
cosmological population models. 

Galactic Halo models 

Galactic Halo models require a halo pop- 
ulation with a large core radius (to avoid an 
anisotropic enhancement towards the galactic 
center) . This is a new population of astronom- 
ical objects, which was not seen elsewhere [p^] . 
It probably require a different distribution (in 
space) than the dark halo material (the lat- 
ter being too concentrated towards the galac- 
tic center). By now there have been several 
suggestions how to form a neutron population 
of this kind. These include either ejection from 
the galactic disk or formation in site. However, 
these models face additional difficulties. 

Approximately 10'^^ ergs are needed for 
bursts at the halo, quite a large amount for 
a neutron star. With a typical size of lO^cm 
this leads to an optical thickness of ~ 10^ 
for 77 — >■ e+e^. These constraints have two 
far reaching implications: First, the optically 
thin neutron star models suggested for galac- 
tic disk sources are inapplicable to grbs at the 
halo. Second, galactic halo sources inevitably 
involve an opaque pair plasma fireball, just 
like cosmological sources ||32|. The physical 
conditions in these fireballs are, however, less 
favorable than the conditions at cosmological 
fireballs for production of grbs. 

Local Population 
Typical objects in the solar system have a 



very small binding energy per baryon and it is 
difficult to imagine a mechanism in which such 
objects generate energies in the 7-ray range 
(see however ||3^). The only hope is probably 
via a magnetic phenomenon. Solar fiares do 
generate grbs which are detected by Compton- 
GRO (these are identified by their location 



and spectrum ^5). However, comparison of 
the size and masses involved in these events 
make it inconceivable that similar conditions 
can be achieved elsewhere in the vicinity of the 
solar system, without leaving any other trace. 

Cosmological Population 



Several groups have shown 

that a cosmological population is compatible 
with the observed V/Vmax distribution. The 
apparent concentration towards us is an arti- 
fact of a combination of redshift effects and 
a possible cosmological evolution. Unfortu- 
nately it is impossible to separate the two ef- 
fects and to determine the typical red shift, 
(Z), to the sources from the V/Vmax distri- 
bution. Depending of the cosmological model 
and the source evolution we have 0.3 < Zav < 



m 



3. For Q = 1 and no evolution Zav ~ 1 

The event rate needed to explain the ob- 
servation is in an amazing agreement with 
the rates estimated for neutron star mergers 
TP| , Because of a historical coincidence 



the forth binary pulsar, PSR1534-I-12, which 
played a decisive role in the determination 
of the merger rate |]TU|, |2D|, was discovered 
39| a few month before Compton-GRO was 



launched and the prediction of the neutron 
star merger rate were not influenced by the 
rates required to explain the Compton-GRO 
results. 

Several other cosmological grbs models 
were suggested after Compton-GRO [0, ^ 



41] , Within the cosmological framework, 
the neutron star merger scenario is the most 
conservative one possible. It is the only one 



based on a source population that definitely 
exists. We know its members will merge, we 
can be certain that huge quantities of energy 
will be released in such mergers, and we find 
the merger rate to be comparable to the ob- 
served burst rate. 

Clues Revisited 

Before concluding we turn once more to 
the clues discussed earlier. The optical depth 
problem disappeared in some sense and re- 
mained in another. Relativistic effects, due to 
the expansion of the fireball []T3|, ^, were not 
taken into account in the original argument 
which is flawed. The resulting spectrum 
from the expanding fireball has the right en- 
ergy range but to a first approximation it is 
thermal. It is a non-trivial (but not impos- 
sible) task to obtain a nonthermal spectrum. 
This problem is shared by all cosmological and 
galactic halo models. 

The March 5th event was one of three soft 
7-ray repeaters, which have a softer spectrum 
and produce repeated bursts from the same 
source, unlike all other sources . It is by now 
generally accepted that these are most likely a 
different phenomenon. 

The nature of the cyclotron lines has 
been fairly controversial since they were first 
reported |§. Mazets etal. g] claimed 
that single "cyclotron absorption lines" were 
present in 20 bursts, with a broad distribu- 
tion of line energies (27-70 keV), but with 
only five lines having energies under 50 keV. 
This is in conflict with the GINGA experi- 
ment which discovered three systems of lines, 
all with nearly identical energies, all under 50 
keV IQ. So far, no lines have been detected 
with any experiment on the Compton Gamma 
Ray Observatory. 

Epilogue: 7-ray bursts circa 2000 



A general test for all cosmological models is 
the expected positive correlation between the 
faintness of a burst (correlated with distance) 
and redshift signatures through the burst du- 
ration and spectrum This correlation 
could be masked by large intrinsic variations 
among bursts, but should eventually be ob- 
served when enough data accumulate. 

At present there are no known optical 
counterparts to grbs. Since neutron star bina- 
ries might be ejected from dwarf galaxies, we 
predict [^], that grbs occur within a few tens 
of arcsecond from dwarf galaxies and within 
but not necessarily at the center of ellipticals. 
Optical identification of some parent galaxies, 
could support this model and the location of 
the burst relative to the galaxy could distin- 
guish this model from other cosmological sce- 
narios that involve supermassive black holes or 
other objects located in the centers of galaxies 
11,11,11. 

The scenario makes one unique prediction: 
strong 7-ray bursts should be accompanied by 
a gravitational wave signal p8|, ^ ^ (though 
the reverse need not necessarily be true if the 
7-rays are beamed). These signals should be 
detected by LIGO when it becomes oper- 
ational (hopefully by the year 2000). LIGO 
should provide good distance estimates to in- 



dividual bursts |43] and should also pinpoint 
the exact time of the merger, in addition to 
furnishing an ultimate proof of this model 
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